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Abstract 

Assuming a double-well bare potential for a self-interacting scalar field, with the Higgs 
vacuum expectation value, it is shown that non-perturbative quantum corrections naturally 
lead to ultra-light particles of mass ~ 10“^^eV, if these non-perturbative effects occur at a 
time consistent with the Electroweak phase transition. This mechanism could be relevant 
in the context of Bose Einstein Condensate studies for the description of cold Dark Matter. 
Given the numerical consistency with the Electroweak transition, an interaction potential 
for Higgs and Dark Matter helds is proposed, where spontaneous symmetry breaking for the 
Higgs held leads to the generation of ultra-light particles, in addition to the usual Higgs 
mechanism. This model also naturally leads to extremely weak interactions between the 
Higgs and Dark Matter particles. 


1 Introduction 

This article aims at providing a dynamical mechanism from which ultra-light scalar par¬ 
ticles arise naturally, and which therefore could be relevant in the context of Bose-Einstein 
Condensate (BEC) Dark Matter (DM), proposed initially in [1] and reviewed in [2]. Axions 
have also been proposed as ultra-light particles to contribute to DM, for which a recent study 
can be found in |3] and a recent review in |Tj. 

Models based on BEC DM assume the existence of a scalar particle for DM, light enough 
for its Compton wave length to be of the order of the size of a DM halo. As a consequence, 
these scalars are in a coherent state, and can be described by a BEC wave function. It is 
usually assumed that condensed particles are non-relativistic, in order to describe cold DM, 
and that they decouple from the Standard Model (SM) at some point in the Early Universe. 

The mechanism proposed here, based on a non-perturbative effect in quantum field theory, 
explains how such ultra-light particle (m <<eV) can arise in the dressed theory, starting from 
a bare theory which contains a typical SM mass scale v. Such an effect is possible if the bare 
theory has several degenerate minima, as the usual double-well potential for a single scalar 
held: since the dressed potential is necessarily convex [5] , quantum corrections must be strong 
enough to eliminate the concave part of the bare potential, which allows for an ultra-small 
ratio m/v. 
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Exploring the possible relevance of this mechanism to BEC DM, we hnd that the pic¬ 
ture is consistent with non-perturbative quantum effects occurring at the Electroweak phase 
transition, where these ultra-light scalars could appear. The stability of such a BEC halo is 
discussed, and the repulsive self-interactions predicted by this mechanism, although tiny, are 
enough to reduce substantially the range of fluctuation wave vectors which could potentially 
lead to a collapse of the condensate, under gravitational forces. Finally, we propose common 
dynamics between the Higgs and the DM particles, which explain how spontaneous symmetry 
breaking for the Higgs held could imply the generation of ultra-light particles, in addition to 
the usual Higgs mechanism. The corresponding bare potential for Higgs/DM naturally leads 
to extremely weak interactions between the Higgs and DM helds, after the non-perturbative 
effect for DM has taken place. 

2 Causality and ultra-light particles 

We review here the features of the convex effective potential (= momentum-independent 
part of the one-particle irreducible generating functional), for a hnite volume. This 4- 
dimensional volume is interpreted as the volume corresponding to the particle horizon, in 
a cosmological context, after requiring that quantum corrections should affect causally re¬ 
lated events only. 

2.1 Convexity of the effective potential 

The convexity of the effective potential for a scalar theory has been known for a long 
time, and is a consequence of its dehnition in terms of a Legendre transform |6]. In the 
situation where the bare potential has several degenerate minima, convexity is achieved non- 
perturbatively, and cannot be obtained by a naive loop expansion. The effective potential 
becomes convex between the two minima of the bare potential, as a result of the contribution 
of several non-trivial saddle points in the partition function [7]. The construction of the 
convex effective potential has been shown explicitly in [8]. We review here the results, but 
the derivation is generalised in the Appendix of the present article, where the real scalar held 
is coupled to a complex scalar held, in order to justify the results presented in section 4. 

We start from the generic double-well bare potential 

Ubare{^) = , ( 1 ) 

and dehne the partition function Z for a hnite space time volume V. The semi-classical 
calculation is done by taking into account both minima of the bare potential, in a saddle 
point approximation, for the calculation of Z. The dressed potential is then, for \(p\ <v and 
in the limit of large space time volume Vv^ » 1, 

Udressedi^) = ^ (^) + ^ (^) + W\ < ^ ■ ( 2 ) 

As expected, this potential is convex, and higher orders in (p are also suppressed by the volume 
(see Figure 1). Outside the minima of the bare potential |(p| > n, quantum corrections are 
perturbative, provided |(^| is not too close to the minimum n, for the “inside” and “outside” 
potentials to match. Note that the dressed potential ([^ is universal in the sense that it 
depends on the bare vacuum expectation value (vev) v only, and not on the bare coupling A. 
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Figure 1: Convexity of the effective potential as obtained in |H], based on a semi-classical 
approximation to calculate the partition function. The effective potential (continuous line) 
is non-perturbative for << where the features of the bare potential (dotted line) are 
suppressed by the volume, as a consequence of the contribution of the two bare vacua. For 
>> only one bare vacuum dominates the partition function, and quantum corrections 
are perturbative. The region ~ is sketched here to match the two previous regimes. 


It has also been shown in |8] that the generalization to an 0(iV)-symmetric scalar theory 
gives, for Vv^ » 1, 


rj(^) 

dressed 


{ 0 ) 


IL(p_ 

2V [v 


7V2 

^ 4(N + 2)V 



+ O(p^) 


for IpI < n , 


(3) 


where p = \/Lp ■ and the case of a complex scalar held is obtained for N = 2. 

The non-perturbative mechanism described here can be understood with tunnelling ef¬ 
fect}^ For a hnite volume, tunnelling between different vacua is allowed, such that the 
partition function is dominated by several saddle points, which leads to a convex effective 
potential. As a consequence, for hnite volume, no spontaneous symmetry breaking occurs. 
It is only in the limit where the volume goes to inhnity that spontaneous symmetry breaking 
can take place, and the system chooses one vacuum among a continuous set of degenerate 
vacua. 

Finally, a comment on the Higgs mechanism might be appropriate at this point. From 
the construction shown in [8] , one can understand that a convex dressed potential cannot be 
obtained if the scalar held if coupled to a gauge held. Indeed, in this situation, one needs to 
hx a gauge in order to dehne the path integral, such that the minima of the bare potential are 
not equivalent, and quantum huctuations are built up from one scalar vacuum only, dehned 
by the choice of gauge. 

In what follows we will focus on the case = 1, since we are looking for a neutral scalar 
held to represent DM. 

would like to thank Arttu Rajantie for this comment 
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2.2 Dynamical generation of ultra-light particles 

In the limit of infinite volume, the potential Q becomes flat between the minima of the 
bare potential. This corresponds to the “Maxwell construction” in Statistical Mechanics, 
arising when two phases coexist in a system. In quantum held theory, this battening means 
that the true vacuum is a superposition of the different bare vacua [5]. This effect has 
motivated studies of inhation [lO], but in this note we will not take the limit of inhnite 
volume. Instead we consider the hnite space time volume of causally related events V = /q, 
where Iq is the particle horizon dehned at a specihc time to to be determined below, and where 
quantum corrections drive the bare potential Q to its dressed form (|^. One identihes then 
the dressed potential ([^ with the expression 


Udressedi,'^'} — 2 ^ ^ 24 ^ ' 

to hnd that the mass m and the coupling constant g are 

1 2 


m = 




9 = 




(4) 

(5) 


From the previous expression, the large volume condition leads to the inequalities 


- « lo « — , 
V m 


( 6 ) 


showing that, necessarily, the dressed mass m is small compared to the bare vev v. We now 
consider the value m ~ 10“^^eV, which is typical in the context of BEG DM mm- Also, 
in order to introduce a typical SM mass scale, we set v to the Higgs vev v = 246 GeV. We 
hnd then that the particle horizon is 


Iq ~ 12 cm , 

(7) 

corresponding to the cosmological time 

to = l/2c^2x 10“^° s , 

(8) 


for a radiation-dominated Universe. It is interesting to see that the choice of the Higgs vev 
for V leads to a cosmological time consistent with the Electroweak phase transition, which 
suggests that the value m ~ 10“^^eV could indeed be related to this transition. Finally, the 
large volume assumption is satished, since vIq ~ 10^^, which leads to 

g ~ 10“®® . (9) 

Although this coupling is extremely small, it contributes in a non-trivial way to the stability 
of the condensate, as explained below, because of the ultra-light scalar mass. We hnally note 
here that detailed models for self-interacting BEG DM are studied in [T3] . 

3 Stability features of the condensate 

We present here two arguments towards the stability of the BEG DM halo, one related 
to the static halo, and the other related to its cosmological evolution. 
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3.1 Gravitational collapse versus repulsive self interaction 

The Authors of HI consider small fluctuations of the non-relativistic condensate, af¬ 
ter linearising the corresponding equation of motion and studying the evolution in time of 
these fluctuations. They conclude that any repulsive interaction between scalars stabilises 
the condensate, but any attractive interaction induces instabilities though low wave vector 
fluctuations. In our present case, two effects compete: the attractive self-gravitation of the 
condensate and the dynamical repulsive self-interaction controlled by 5 ^ > 0 . 

It is explained in [H] that, if one takes into account gravity only, fluctuations with wave 
vector k evolve in time as exp(± 7 o(/c)f), with 


k IQiiGnm? 


( 10 ) 


where n is the number density of scalars of mass m and G is the Newton constant. It is easy 
to see that 70 vanishes for the critical wave vector k^ given by 




70 

m 


m 




( 11 ) 


• For k > /cq, 7 o is purely imaginary and fluctuations are purely oscillatory: the homogeneous 
conhguration of the condensate is stable; 

• For fc < fcg, 7 o is real and the condensate is destabilised by low fc-modes: the condensate 
collapses. 

The typical mass density of our galaxy is 


Pgai ~ 1 GeV cm ^ 

and the number density of scalars is thus 

1 GeV cm“^ 


n ~ 


~ 10^° eV^ , 


10-23 eV 

such that, together with G ~ 10“^® eV- 2 , one hnds 

h* 

^ ~ 10-3 , 
m 


( 12 ) 


(13) 


(14) 


The fluctuation wave lengths ~ k~^ which could destabilise the condensate are therefore much 
larger than the BEG halo radius ~ m~^. If one takes into account the repulsive interaction 
controlled by the coupling constant ([^, 70 (/c) is replaced by 


,,, k / on IGirGnm^ 


(15) 


and the critical wave vector k* which satishes 7 (A:*) = 0 is now given by 


gn 

4m3 


—1 -|- W1 -|- 


2567rGm^ 


g'^n I 


(16) 
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Numerically, one finds 


k* 

— ~ 10 
m 


-16 


h* 

ISjn 

« — 
m 


(17) 


which reduces enormously the range of modes which potentially could destabilize the con¬ 
densate. A more detailed analysis would be required, going beyond the linearised equation 
of motion for fluctuations, to conclude on the definite stability of the condensate. But one 
can see here that the repulsive interaction predicted by the present model provides a huge 
improvement to the BEC halo stability. 


3.2 Time evolution 

The dressed potential (|^ holds as long as |(p| < n, but this condition remains valid for 
later times t > to, since the energy density decreases. More precisely, the equation of motion 
for the field is 

(/3 -f 3Hip -|- = 0 , (18) 

where the Hubble parameter is if = for the radiation-dominated era. We neglect here 

the back-reaction of the scalar’s dynamics on gravity, which is assumed to be determined 
by radiation. Since if (to) — 4 x 10“®eV>> m, one can initially neglect the mass term in 


eq.(18), which can be integrated once to give 


lf = lfO 


3/2 


(19) 


where ipo is a constant, corresponding to an initial condition. The energy density is then, if 
one neglects the potential term /2, 


P = 




( 20 ) 


and is proportional to the 6th inverse power of the scale factor a{t) = consistently with 
The approximation (20) is not valid anymore around the time ti, when the Hubble 


parameter H{ti) becomes equal to m. This time is 


~ 3 X 10^ 


( 21 ) 


which is still during the radiation-dominated era. A more detailed analysis can be found 
in [12], showing that the energy density goes on decreasing, such that the amplitude of the 
field remains smaller than the bare vev v, and the potential ([^ remains appropriate for the 
description of non-interacting ultra-light particles. 


4 Higgs/DM dynamics 

We suggest here common dynamics to the Higgs and the ultra-light DM particles, given 
the above numerical coincidence. The motivation is to generate ultra-light particles simulta¬ 
neously with the Higgs mechanism. As shown here, the non-vanishing vev of the Higgs field 
could also generate the double-well bare potential for the DM neutral field, which therefore 
would lead to ultra-light scalars, after non-perturbative quantum corrections. 


6 




For simplicity, we consider the Abelian Higgs model, and extend it with the nentral DM 
scalar tp 

c = + D^<p{Dyy + - Ut,are{<P<P*,p") , ( 22 ) 


where 0 is the complex Higgs field and is the covariant derivative. The potential is chosen 
as 


Uhare 






( 23 ) 


which is minimised for the field values 




( 24 ) 


The derivation of the semi-classical effective potential for DM is given in the Appendix, which 
generalises the original derivation in [8] to the situation where the real scalar field is coupled 
to a complex scalar field. The result found in the Appendix can be obtained in a simpler way 
though, as follows. Because the Higgs field is coupled to the gauge field, its vev v is unique 
and is determined by spontaneous symmetry breaking. But this is no the case for the real 
field p which, for a given Higgs field configuration, sees the two bare minima Fy/cjxjF. The 
non-perturbative mechanism described in section 2 leads then to the DM effective potential 


U^ff _ 
DM ~ 


1 




+ 




2V 00 * 121 / (00 




+ 0{py for p^ < 


( 25 ) 


We then parametrise the Higgs field in the usual way 




( 26 ) 


where ^ is the dynamical Higgs field and the Goldstone mode rj is absorbed as a longitudinal 
component of the vector field, and is therefore ignored in what follows. An expansion of 00* 
around to the second order in y leads to 



2V [vJ 12V [vJ 


V v3 2V 


( 27 ) 


where dots represents higher-order interactions terms. This result is consistent with the full 
derivation given in the Appendix, and describes the following Higgs-DM interactions: 


• A decay process Higgs—?-DM-DM which, given the ultra-small ratio m/mniggs — 10“^^, 
leads to relativistic DM particles, not contributing to the BEG. These relativistic DM 
particles are predicted in BEG DM models [2], and consist in a small portion of DM 
particles. This cubic interaction also represents an annihilation of two DM particles 
into a Higgs particle, which is kinematically possible with relativistic DM only; 




A scattering process Higgs-Higgs—)-DM-DM or Higgs-DM—)-Higgs-DM, which is repul¬ 
sive and therefore stabilises further the DM halo, according to the analysis of section 


3.1. 


These two interactions are highly suppressed though, by coupling constants which are pro¬ 
portional to (l/n^)“^ ~ 10“®®, and therefore are not detectable. 

This conclusion holds as long as the amplitude of p"^ is smaller than where DM particles 
see the effective potential (27). In a regime where the number density of DM particles is 
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such that (p^ > though, the potential seen by these particles consist in perturbative cor¬ 
rections to the bare potential (23), such that the interaction DM/Higgs could be detected. 
The critical held value characterising the transition between perturbative and non- 

perturbative effective potentials (see Figure 1) corresponds to an energy density of the order 
0) ~ V~^, and therefore to the number density of DM particles 


^crit 


_ n _ 246 GeV 
m /q (12 cm)2 


~ 10*^^ eV' 


The ratio between ricru and the number density (13) of DM in the galaxy is huge 


"^crit _ ^ q 37 
"^gal 


(28) 


(29) 


but can be achieved during the quark-gluon plasma phase generated in heavy-ion collisions. 
Indeed, the critical energy density for a quark-gluon plasma to be created corresponds to the 
proton energy density 


Pproi - 1 GeV fm ^ 


(30) 


whose ratio with the galaxy energy density (12) is 


^ . (31) 

Pgal 


An experimental signature for the present model through the study of quark-gluon plasma 
might therefore be possible, although much more work needs to be done in this direction. 


5 Conclusion 

The main point in this article consists in the identifications (|^, which allow the hierarchy 
between the scales m and v in natural way, and predict that ultra-light particles have repulsive 
self-interactions which improves the stability of the BEG DM. The consistency of the picture 
with the Electroweak transition shows a new avenue to explore: the link between the Higgs 
particle and Dark Matter. Both particles could see the same vev, but behave in a different 
way, according to their coupling to gauge fields: the Higgs field sees perturbative quantum 
corrections, whereas the additional scalar, blind to gauge fields, undergoes non-perturbative 
corrections and becomes ultra-light. The common dynamics between the Higgs and the DM 
fields, presented in section 4, is a first step in the direction of a potential unification of the 
two fundamental origins of mass in the Universe. 


Appendix: Semi-classical derivation of the DM effective 
potential 


We generalise here the derivation given in [S] to the situation where the real scalar field 
is coupled to a complex scalar field which sees only one bare vacuum because of gauge fixing. 
Starting with the Abelian model (22), we focus on the scalar sector. As shown below, the 
main features of the mechanism are: 










• Because of gauge fixing, the Higgs field sees only one vacuum, and the corresponding 
semi-classical approximation doesn’t modify the bare Higgs sector; 

• The real scalar field sees two bare vacua, since it is not coupled to the gauge field. 
Tunnelling effect - for finite volume - leads then to the non-perturbative mechanism 
which results in the suppression of the effective potential by the volume. 


We follow here the usual steps of path integral quantisation, where the Euclidean partition 
function, functional of the sources j,j*, k, is 


Zli, f, ^] = I V] exp j + Ui^re{p, v)+]<!> + fV + kip'j , 

(32) 


where p = and 


Ubareip, = || - v'^Y + {p^ - 


2\2 


(33) 


We choose the (unique) Higgs vacuum to be in the real direction and disregard the Goldstone 
mode, such that 0* = 0 and the partition function reads 


Z[r,k] = / T>[0,(p]exp 


ip + Ubareip, + ‘2rp + kip] ) (34) 


where r = \j\. Due to the double-well structure of the potential, this partition function is 
dominated by two uniform saddle point configuration (p+, 0+) and (p_,0_). These saddle 
points are obtained after minimizing the action, for uniform sources, and are solution of the 
equations 

^oipl-pl) + k = 0 and poi2pl - pi - 1) + r = 0 , (35) 

where the dimensionless helds are 


Po 


V 


Po 

, Po = — 

V 


12 r 

r =- 

9hv^ 


~ 6k 

k = -: 

9hv 


(36) 


and Po = p±, Po = P±- A Taylor expansion gives for the hrst saddle point 


= i-r/2- k/2 - 3r^/8 - 3fk/A - 5p/8 - /2 - 3f‘^k/2 - 9hP/4 - 25p/16 

-105hVl28 - 105h^fc/32 - - 69hP/8 - 637fcVl28 + ■ ■ ■ 

p+ = l-f/2-k- 3hV8 - 5hfc/4 - 33p/3 - 7^2 - 9h2fc/4 - 75hP/16 - 4p 

-1057V128 - 1477^^/32 - 2077^P/16 - 6377^/32 - lOSfcVS + ''' (37) 


and for the second 

p_ = 1 - 7/2+ fc/2-37^/8+ 37fc/4-5p/8-7V2 + 372 ^/ 2 -97p/4 + 25p/16 
-1057V128 + 1057^fc/32 - 4417^^/64 + 697p/8 - 637fcVl28 + ■ ■ ■ 
p_ = -1 +7/2 -fc + 37V8-57fc/4 + 33p/3 +7^2- 972^/4+ 757 P/I 6-4P 

+ 1057V128 - 1477^fc/32 + 2077^^/16 - 6377p/32 + lOS^VS + ■ ■ ■ (38) 

^This is true as long as the sources are smaller than critical sources, or equivalently the classical fields 
(defined below) are smaller than the bare vacuum v, which is what we assume here. As explained in [5], 
for classical fields larger than v (or equivalently for large sources), the partition function is dominated by 
one configuration only, and the semi-classical approximation that we use below, to determine the effective 
potential, simply leads to the bare potential for these field values. 
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where dots represent higher orders in the sonrces. The partition fnnction is then determined 
by the saddle point approximation 

Z[r,k\ ~ ^exp{V[Ubare{p+,^+)+ ‘^rp++ k(f+]) 

+ ^exp(y[Ubare{p-,P>-)+ ‘2rp_ + kif-]) , (39) 

where V is the space time volnme. A Taylor expansion in the sonrces gives 

Z[f,k\ = 1 + 4Ah + 8AV + + (32AV3)h^ + 32A3hP 

+ (32AV3)h^ + + (32AV3)P + ■ • • (40) 

where 

, _ gnVv^ 

24 

and only the dominant terms in A >> 1 are kept for the large volnme approximation. For 
a finite volnme and constant field confignrations, the fnnctional derivatives become partial 
derivatives and the classical fields are given by 



1 


1 

dZ 


5r 

^ iAZ 

dr 


1 6Z 

^c = v(pc = ^ ^c = 


-1 dZ 


4AZ dk 

In the large volnme approximation A >> 1 we find the Taylor expansions 


pc = 1 - (l/2)h-(3/8)P-2AP - (1/2)P-2AhP 
-(105/128)P - (15A/4)Pp - (32AV3)P + ■ ■ ■ 
pc = —4Afc + 4Arfc + 2Ar^fc + (64A^/3)P 

+ (5A/2)Pfc + (128AV3)hP + • • • , 


(42) 


(43) 


where one can see that, whereas the complex classical amplitnde pc oscillates aronnd the vev 
n, the real scalar pc oscillates aronnd 0 becanse of tnnnelling effect. The latter relations can 
be inverted to obtain the sonrces as fnnctions of the flnctnations = P — 1 and pp- 

r = -2|e-3el-^p(4A)-a + 3e>P(4A) 

-3ecVP(2A)-^P(12A) + --- 
k = Pc/{AA)-lPc/{2A) + ?>ilPc/{AA)+pl/{l2A) 

-llVc/A - lcPll(P>A) + ■■■ (44) 


where only the dominant terms in A >> 1 are considered. Note that, in the present semi- 
classical approximation, the flnctnations are negative (which can be see from eq.(43)), such 
that r in eq.(44) is positive, as expected, since it is the modulus of the source j. 

The effective action for uniform fields 


V[ic,ifc] = VU^f^{ic.Vc) , (45) 

is defined as the Legendre transform of the connected graphs generating functional 

W[r,k] =-\xiZ[r,k] , (46) 


10 













and reads, for uniform classical fields, 


r[ee, Vc] = W^[r, k] - Vi2ri, + kif,) , (47) 

where the sources r, k have to be understood as functions of the classical helds through 
the relations (44). From its dehnition, it is known that F satishes the two following equations 


1 

6T 

1 

dr 

r = vr = - 
2 


AA 

die 


k = vk = 


6r 

6ipc 


k = 


-1 dr 

AA d(pc 


which, together with the expansions (44), can be integrated to get 

F = AAll + AAll + All + + \pI^I + 

9hVv^ 


24 


(1 + ^cY 


I ~2 I ~4 ~2 I ~2 ~2 I 

+ 2^c + ~ + 2dc7^c + 


The effective potential is finally obtained after dividing by the volume V 




and contains two contributions: 


1 

2 \v J ^12Vn 




+ X 


2 


+ 


(48) 


(49) 


(50) 


• The part proportional to gjf, which is the usual bare Higgs potential, and has not been 
modihed in this semi-classical approximation; 


• The part depending on obtained in eq.(27) with a simpler argument, which is sup¬ 
pressed by the volume. 
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